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In a previous paperEl a method was developed to determine the shapes of 
dark matter halos of spiral galaxies from the flaring and velocity dispersion 
of the gas L awn Here I present the results for the almost edge-on Scd galaxy 
NGC 4244E3E3 and preliminary results for the Milky Way. NGC 4244's 
dark matter halo is found to be highly flattened with a shortest-to-longest 
axis ratio of 0.2^' j. If the dark matter is disk-like, the vertical velocity 
dispersion of the dark matter must be ~ 20% larger than the measured 
tangential dispersion in the H I . The flaring of the Milky Way's gas layer, 
the local column of identified stars and the total column within 1.1 kpc from 
the plane are consistent with a moderately flattened dark halo (E7-E0) and 
galactic constants of Ro = 7.1 kpc, Go = 180 km s _1 , while the stellar disk 
is ~ 80% of maximal. 



1 Introduction 

Although rotation curves of spiral galaxies have been used as evidence 
for the presence of dark matter (DM), little is known about the nature, 
extent and actual distribution of the DM in individual galaxiesoEj. The 
equatorial rotation curve, which probes the potential in only one direc- 
tion, provides no information about the shape of DM halos. 

Several methods have been usedto determine the shapes of dark 
matter halos. The analysis of warpal^ shows that only one of the five 
systems studied requires a DM halo as flattened as E4||. On the other 
hand, in studies of polar ring galaxiesc§E3Ej substantially flattened DM 
halos are found (E6-E7 and E5). The shape of the dark halo of the Milky 



a To appear in the proceedings of "Aspects of Dark Matter in Astro- and Par- 
ticle Physics", Heidelberg, September 1996 (World Scientific; Eds. H.V. Klapdor- 
Kleingrothaus, Y. Ramachers). 

b A spheroidal system with shortest-to-longest axis ratio c/a (q p ) of the density 
contours has a shape En, with n such that q p = 1 — n/10 
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Way has been estimated from stellar kinematics (EO - E7)oOE3H, and 
the method of flaring gas layers (E0-E7; see below). The SO galaxy NGC 
4753 seems to have a rather round DM haloi (El). The X-ray isophotes 
of two early type galaxies indicate that their dark halos are moderately 
flattened! (E5.5 and E6). Cold dark matter galaxy formation simula- 
tions which include gas dynamics tend to produce rather oblate DM 
halosll, with q p = c/a = 0.5 ± 0.15 (Fig. 1). 

Fig. 1. A histogram 
of the known DM halo 
shapes. The Dot- 
ted line and the filled 
squares represent the 
theoretical prediction. 
The points with error 
bars represent the in- 
dividual galaxies. 
Note the discrepancy 
between the results 
from the warping-gas- 
layer method (right- 
most bin) and the 
other methods. 




0.4 0.6 
DM axis ratio c/a 

The shape of the dark halo can be determined by comparing the 
measured thickness of the gas layer (flaring) with that expected from 
a self-gravitating gaseoiisdisk in the potential due to the stars and 
the flattened DM halcc3E3. If a round halo, with a certain density 
distribution, is squeezed along the vertical axis, the densities and the 
exerted gravitational forces will increase, resulting in a thinner H I disk 
and higher rotation speeds. In order to fit the observed rotation curve, 
one has to deform the DM-halo such that the DM-halo density (pdm) 
at large distances will be roughly inversely proportional to the halo 
flattening q p . Since the thickness of the gas layer beyond the optical 
disk is proportional to 1/y/pDM, Q'pOC (width of the gas layer) 2 . 

Below I apply the method to the galaxyNGC 4244 for which the 
basic parameters were determined in Paper IlE3, and present preliminary 
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results for the Milky Way. 



2 The Method 

For a constant temperature gas, the density at height z above the plane, 
Pgas(z), can be calculated from the equation of hydrostatic equilibrium: 

with Ugas the velocity dispersion. Assuming that the density distribution 
p(z) extends to infinity (i.e., p(R,6,z) = p(z)) the potential has only 
a vertical gradient so that the vertical force can be calculated from the 
Poisson equation: K z (z) = —4ttG Jq p(z')dz' . Although not perfect, 
this plane parallel sheet approximation has been used extensively in the 
past and can be used faLfirst order approximations to the expected 
thickness of the gas layered. For example, in the cases that the gaseous 
self-gravity, the stellar disk, or the DM halo dominates the potential the 
full width at half maximum (FWHM) of the gas layer is given by: 

FWHMgaJ (R) « 1.58 a gas /H gas (2) 

FWHM S g ir(R) ~ O.GyJzjEstars X O gas (3) 



FWHM^°(R) « 2.35, 



2Aq p \ fa 



\ \1A + q p ) U 



'gas 



R c {q p f + R? (4) 



with Y, gas and X s tars the gaseous and stellar surface densities in units 
of M Q pc~ 2 , z e the scale-height of the stellar disk, V max the maximum 
rotation speed (in km s" 1 ), R c (q p ) the core radius of the dark halo, and 
R the cylindrical radial distance (all distances in kpc). These equations 
can be combined to yield an approximate value for FWHM gas if more 
than one component contributes to the potentialc3. In regions where 
the rotation curve rises steeply or the surface density varies rapidly this 
approximation fails so that it is better to calculate the vertical force 
from the density distribution of the whole galaxyc3: 

roo roo fir J 

K z (R,z) = G rdrp{r,0) dwp tot {r,w) — — -=-, (5) 
Jq J-oo J-7T dz s — o 
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with s = {r, 8, w}, S = {R, 0,z}. I incorporate three components in the 
global mass model; 1) a stellar disk with a density distribution which 
decreases exponentially with radius and height above the plane; 2) a 
non-singular flattened isothermal DM-halo with core radius R c (q p ) and 
central density po(q p ) with a density distribution proportional to 1/R 2 
to reproduce the observed "flat" rotation curves; and 3) a gaseous disk. 

Of course, the true DM-halo density distribution may be different]!!. 
However, for roundish DM distributions the vertical force is roughly pro- 
portional to the radial force which is the same for all disk-halo combina- 
tions that reproduce the observed rotation curve, so that at large radii 
the width of the gas layer is independent of the radial distribution of 
the DM. The flattening of the DM halo introduces a ~ y^-dependence 
on the thickness of the gas layer. 

Comparing the thickness of the gas layer beyond the optical disk 
with model flaring curves, calculated for a series of models with varying 
halo flattening, then yields the halo shape. 



3 Results: NGC 4244 

The Scd galaxy NGC 4244 was observed for 14 hours with the VLAj] to 
determine the gaseous velocity dispersion and the rotation curve. While 
compact —fast rotating galaxies are known to have declining rotation 
curvest9'c3, NGC 4244 is the only low mass galaxy (V max « 100 km s _1 ) 
for which the rotation curve falls, in Keplerian fashion. 

A new technique^ has been used to determine the width of the gas 
layer (Fig. 2). An upper limit to FW H M gas is found by assuming a 
constant inclination of 84?5, while incorporating the slight warp into the 
analysis yields the best values for the flaring (filled triangles). 

Comparing the observations with the model flaring curves (drawn 
lines), I conclude that the DM halo of NGC 4244 is highly flattenedH 
q p = 0.2 ± 0.1 . This flattening, the most extreme value reported to 
date, lies at the extreme end of the theoretical predictions (Fig. 1). Is 

c The VLA of the National Radio Astronomy Observatory is a facility of the Na- 
tional Science Foundation operated under cooperative agreement by Associated Uni- 
versities, Inc. 
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it possible that some systematic effect plays a role and that NGC 4244's 
DM halo is less flattened? I investigate several possibilities: a) The 
systematics introduced by the uncertainty in the inclination is discussed 
in the caption of Fig. 1. b) Non-thermal pressure gradients due to 
magnetic fields and cosmic ray heating [which add to the LHS. of Eqn 
(1)] are unlikely to be important beyond the optical disk because cosmic 
rays are closely related to sites of star formations. If they are 

Fig. 2. The measured 
gas layer widths (open 
and filled triangles for the ~ , 5 
i=84°5 and warp-included - 
case, respectively). The | 1 
drawn lines correspond to 05 
models with different halo 
shape:g p =1.0, 0.7, 0.5, 0.3, 
0.2, and 0.1, from top to 2 
bottom. Using both inch- ~ 

o 1.5 

nation cases I find: q p = i 
0.2±g;f . The mass-to-light i i 
ratio of the stellar disk is 
not constrained by these 
flaring measurements. 

° 5 10 

Galactocentric radius (kpc) 

important an even denser, i.e. flatter, DM halo would be required to 
have a gas layer as thin as observed. And c) If the vertical velocity 
dispersion is smaller than the measured planar dispersion, the DM halo 
would be rounder than inferred above. There is no observational evi- 
dence that such might be the case. Furthermore, because the interstellar 
medium is likely to be in the warm neutral phase due to the low pres- 
sure^, the short collision times (< 10 5 year) preclude any anisotropy in 
the velocity dispersion tensor. We conclude that the DM halo of NGC 
4244 is significantly flattened, with q p — 0.2^^ . 
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3.1 An Alternative Explanation ? 



Pfenniger et aZ.H reviewed disk-like molecular hydrogen as a dark mat- 
ter candidate. The fact that in many galaxies the shape of therotation 
curve due to the gas is similar to the observed rotation curvaffl could 
then be explained if ~ 6% of the gaseous surface density is in atomic 
form. Since the self-gravity of the gas layer beyond the optical disk 
strongly affects the flaring (Eqn 2) I investigate whether the cold gas 
hypothesis is consistent with NGC 4244's flaring curve. In order for the 
dark disk to have a thickness equal to the H I layer, the vertical velocity 
dispersion of the dark disk must be ~ 20% larger than the dispersion in 
the H I . Furthermore, to avoid radial instabilities, the planar velocity 
dispersion has to be 20% - 100% larger than the vertical dispersion: a 
dark disk requires an anisotropic velocity dispersion tensora. 



4 Results: The Milky Way 

The flaring of the Milky Way's gas layer can be used to constrain the 
shape of the dark matter halo as well. The analysis is complicated 
however because the uncertainties in the distance to the galactic center 
(Ro ), the rotation speed at the solar circle (Go ), and the scale-length 
of the stellar mass (/ir). Sackett£3 reviewed these values and found 
the following: R = 7.8 ±0.7 kpc, 9 = 200 ± 20 km s _1 , and /i R = 
3 ± 1 kpc. Due to independent constraints on the Oort constants^ 
A — B (= 00 /Rq ) and A, small values for Rq requires small values for 
Go , and vice-versa. I calculated flaring curves for Milky Way models 
with a small, intermediate and large value for Rq : (Rq ,Qq )=(7. 1,180), 
(Rq ,0 o )=(7.8,200), and (Rq ,0 o )=(8 f 5,220). These Milky Way mod- 
els include a stellar bulge and disked which is truncated at around 
(Rq +4.5) kpctil. Various authors re porte d on the thicknes 
and column- and volume densitiesBcHlaEa of the H I & H2 layers of the 
Milky Way. These observations were all scaled using a fit to the rota- 
tion curve of the innere3 and outerca Milky Way, scaled separately for 
the three Rq & ©o combinations. Since the models which predict the 
thickness of the gas layer cannot handle a multi-component interstel- 
lar medium, the Milky Way's flaring data must be compared with the 
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model predictions beyond the solar circle, where the H2 column densi- 
ties are small. In the model calculations I used the velocity dispersions 
of the H I (9.2 km s _1 ) as determined in the inner galaxycEl. 

Fig. 3 compares the observed thickness of the gas layer with model 
predictions for various halo flattenings. Only those model flaring curves 
are plotted for which the column density in identified starsl3 (35 ± 5 
M@pc~ 2 ) and the total column density within 1.1 kpc from the planalZI 
(71 ± 6 M Q pc~ 2 ) are within 3-cr of the nominal values. 




Galactocentric radius (kpc) 

Fig. 3. The flaring of the Milky Way (triangles) superposed 
on model flaring curves with different halo shapes. For this as- 
sumed iio -©o combination, the rotation curve has been deter- 
mined till 18 kpc, and extrapolated beyond. The last three data 
points are thus least reliable. This disk model is 70% (7 = 0.70) 
of maximal. The errors are not random errors, but are calcu- 
lated by averaging the measurements of various authors (3 for 
R < 18 kpc, 2 beyond). The three values plotted to the right of 
the model curves indicate: a) the halo flattening, b) the differ- 
ence between the model and observed (stars+gas+halo) surface 
density density within 1.1 kpc, and c) the difference between 
the model stellar column density and the column of identified 
stars (in units of their rms uncertainties) . 

Model calculations with an optical scale-length of 2 kpc yield a some- 
what better correspondence between model and flaring data; longer 
scale-lengths fit somewhat worse. The same is true for Milky Way mod- 
els with larger values for Rq . Because the errors on the observed widths 
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are systematic, no formal determination of the dark halo flattening can 
be made. For Rq =7.1 kpc, halo flattenings between 0.3 and 1.0 are con- 
sistent with the observed flaring. For larger values of Ro , the correspon- 
dence between the model curves and observations become progressively 
worse, while the inferred value for q p increases with Rq as well: q p = 
0.75 ± 0.25 for R =7.8 kpc. For 9 = 220 km s~\ the model flaring 
curves lie below the observations for all choices of Rq , halo flattening 
and optical scale-length. 

The inferred halo flattening depends only slightly upon the actual 
values of the scale- length of the disk and its M/C because the halo flat- 
tening; influences the thickness of the gas layer only beyond the optical 
disked. While mass models with /ir=3 kpc require stellar disks which 
are (60 ± 10)% of maximal, a disk with an optical scale- length of 2 kpc 
requires 7 ~ 0.9 ± 0.05. 

For any given mass model, the midplane dark matter mass density 
at the solar circle is found to be proportional to gr-°- 72±0 - 02 ) with an 
average value of: p(z = 0; q = 1) = (9.4 ± 2) 10~ 3 M p C - 3 = i0.38 
± 0.08) protons cm~ 3 . Due to the strong (/^-dependence of R c {q p f^, the 
halo midplane density depends more strongly upon q p near the center 
of the Galaxy than at large radii (where pdm{z = 0) oc q~ 5 ). 

This determination of the flattening of the Milky Way's dark halo 
(0.65 ± 0.35) is consistent with previous determinations!! 3 !!!! and in- 
sensitive to the stellar mass distribution. The major uncertainty is the 
value of the gaseous velocity dispersion: a 15% smaller dispersion would 
yield a round halo. These preliminary calculations favour Galaxy mod- 
els with low rotation speed and a small distance to the Galactic center. 

5 Looking Ahead 

I have presented the results of a new method to determine the shape 
of dark matter halos from sensitive H I measurements and careful mod- 
eling. The first results exclude neither cold dark matter nor disk-like, 
baryonic dark matter. With current technology and "reasonable" ob- 
serving times, the thickness of the H I layer can be measured for galaxies 
closer than ~15 Mpc at inclinations ^ 60°. I recently observed seven 
more systems (NGC 2366, 2403, 2903, 2841, 3521, 4236, and 5023) for 
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which I will try to determine the DM halo shapes. Furthermore, the 
analysis of the flaring of the gas layers of M31 is in progress. With 
this increased sample it will be possible to gauge the significance of the 
highly flattened halo of NGC 4244 and, hopefully, put more stringent 
constraints on the nature of the dark matter. 
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